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^ ' Abstract. Numerical homotopy continuation of solutions to polynomial equations is 

the foundation for numerical algebraic geometry, whose development has been driven by 
applications of mathematics. We use numerical homotopy continuation to investigate the 
problem in pure mathematics of determining Galois groups in the Schubert calculus. For 
example, we show by direct computation that the Galois group of the Schubert problem 
of 3-planes in meeting 15 fixed 5-planes non-trivially is the full symmetric group 56006- 
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'nI" ■ solve systems with over 40 million solutions |T9]. The emerging field of numerical algebraic 

^ geometry [251 126] uses numerical homotopy continuation as a foundation for algorithms 
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Introduction 

Numerical homotopy continuation [26j gives a method to find all solutions to a system 
of polynomials with finitely many solutions. Current parallel implementations [18] can 



O . to study algebraic varieties. While numerical algebraic geometry was developed for ap- 

plications of mathematics, we apply it in pure mathematics, computing Galois groups of 
enumerative-geometric problems from the Schubert calculus, called Schubert problems. 
Along with [2] , this is one of the first applications of numerical algebraic geometry to a 
! problem in pure mathematics. 

Jordan introduced Galois groups of enumerative problems in 1870 [11] and Harris laid 
their modern foundations in 1979 [8], showing that the algebraic Galois group is equal 
^ . to a geometric monodromy group. Byrnes [H Section 5] used Harris's theory to prove 

that the general problem of placing poles with static output feedback in linear systems 
theory was not solvable by radicals. He used numerical homotopy continuation to show 
that a particular Galois group arising in pole placement was the full symmetric group, 5*5. 
Underlying this calculation was a Schubert problem. Vakil [52] applied his geometric 
Littlewood-Richardson rule [21] to study Galois groups of Schubert problems and showed 
that many Schubert problems have Galois group containing the alternating group. 

A Schubert problem is simple if it involves no more than two Schubert conditions of 
codimension more than 1. Simple Schubert problems are natural to study [28l 129] and 
among all Schubert problems on a given Grassmannian, they have the largest intersection 
numbers, so they are the most challenging for direct computation. They may also be 
formulated as complete intersections, which is a restriction imposed by our software. 
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Numerical Theorem. The Galois group of the Schubert problem of 3-planes in 
meeting 15 fixed 5-planes non-trivially is the full symmetric group 5*6006 ■ 

This is a numerical theorem, as our software does not certify its output. We have 
computed Galois groups of scores of other simple Schubert problems, including one on 
the Grassmannian of 4-planes in C*^ having 8580 solutions and ones on Grassmannians of 
3-planes in and in having 10329 and 17589 solutions, respectively. In every case, we 
find that the Galois group is the full symmetric group. Table[2]in Section [3^ records some 
of these calculations. Based on this evidence, we conjecture that every simple Schubert 
problem on a Grassmannian has Galois group equal to the full symmetric group. 

Not all Schubert problems have Galois group equal to the full symmetric group. Using 
an idea of Derksen, Vakil [22] gives some Schubert problems on Grassmannians whose 
Galois group is not the full symmetric group and Ruffo, et. al. |2Ij give one on a particu- 
larly small flag manifold. None of these examples can be studied with our software, which 
requires the Schubert problem to have a formulation as a complete intersection. 

Our software has two implementations in Maple which use homotopy continuation to 
compute elements in the Galois groups and either Maple or GAP |6] to determine if these 
elements generate the full symmetric group. For the continuation, both implementations 
use PHCpack |33j through its Maple interface PHCmaple [17j and the second may also 
call Bertini pQ. The advantages of Bertini are that it can use arbitrary precision and it 
gives an independent verification of our results. 

Numerical techniques give insight into some mathematical properties that are far be- 
yond the reach of other methods. For example, Billey and Vakil [3] studied Galois groups 
of Schubert problems using symbolic methods. The largest problem that they treated 
(showing its Galois group is the full symmetric group) had 9 solutions on the Grassman- 
nian of 2-planes in C^, and they stated that the Schubert problem on this Grassmannian 
having 14 solutions was computationally infeasible. 

The largest simple Schubert problem which we have solved symbolically has 91 solu- 
tions pn §5.3]. In contrast, numerical methods allow us to solve Schubert problems with 
as many as 17589 solutions. These examples actually underestimate the gap between 
the computational possibilities of symbolic and numeric methods, because they were per- 
formed on serial machines. 

Current and (likely) future increases in computer power will come from multiple core 
and distributed computing. This is a break with the past, when improvements in compu- 
tational power came from increasing the clock speed of single-processor units. Symbolic 
algorithms have limited potential in this regime, as Grobner basis computation appears to 
be intrinsically serial and thus can not be efficiently parallelized. In contrast, numerical 
homotopy continuation is easily parallelized, since its atomic tasks are independent. Thus 
methods based on numerical continuation will reap the benefits of future parallel archi- 
tectures. In addition, numerical algorithms typically require less memory than symbolic 
algorithms. In particular, the sizes of final and intermediate expressions in Grobner basis 
computation not only may be large, but also are unpredictable. For these reasons, we feel 
that the future of computing in algebraic geometry lies in numerical algorithms. 
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This paper is structured as follows. In Section [T], we describe the basic geometry of 
Schubert problems and Harris's theory of Galois groups. In Section [21 we explain the use 
of homotopy continuation for simple Schubert problems. We present our software and 
algorithms and discuss our results in Section 3, which include the computation described 
in the Numerical Theorem. We describe future work in Section 4. 

1. Galois group computation of Schubert problems 

The Schubert calculus [14j is a method to compute the number of solutions to Schu- 
bert problems, which are a class of geometric problems involving linear subspaces. The 
prototypical Schubert problem is the classical problem of four lines: How many lines in 
space meet four given lines? To answer this, note that three lines ^1,^*2, ^3 lie on a unique 
doubly-ruled hyperboloid, depicted in Figure [1] These three lines lie in one ruling, while 




Figure 1. The two lines meeting four lines in space. 

the second ruling consists of the lines meeting the given three lines. The fourth line £4 
meets the hyperboloid in two points. Through each of these points there is a line in the 
second ruling, and these are the two lines mi and m2 meeting our four given lines. 

The Galois group of this Schubert problem is the group of permutations which are 
obtained by following the solutions over loops in the space of lines ^1, . . . ,^4. Rotating 
£4 about the point p gives a loop which interchanges the two solution lines mi and m2, 
showing that the Galois group is 6*2, the full symmetric group on two letters. 

1.1. Schubert problems in the Grassmannian. A typical Schubert problem asks for 
the linear subspaces of a fixed dimension [k-planes) in C" that have specified positions 
(incidence conditions) with respect to some fixed, but otherwise general, linear subspaces. 
Each incidence condition defines a set of /c-planes, called a Schubert variety, and the 
solutions to the Schubert problem are the points of intersection of the corresponding 
Schubert varieties. We describe this class of problems. 
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The Grassmannian G{k,n) is the set of /c-planes in C". This is a complex manifold of 
dimension k{n—k). The problem of four lines involves the four- dimensional Grassmannian 
G{2, 4) as a line in (projective) 3-space corresponds to a 2-plane in C^. The set of lines m 
meeting a fixed line i corresponds to the set of 2-planes M of C"^ whose intersection with 
a fixed 2-plane L is at least one-dimensional, and this set of lines m is a Schubert variety. 
The problem of four lines asks for the points common to four such Schubert varieties, one 
for each of the given lines ii — £4 in projective 3-space. 

The specified positions of /c-planes in Schubert problems are in reference to flags in C". 
A flag F, is a sequence of linear subspaces 

F. : Fi C F2 C ■ ■ ■ C C F„ = C" , 

where i = dimFj. The possible positions are encoded by partitions. A partition A is a 
weakly decreasing sequence of integers 

A : (n-k) > Ai > A2 > ■ ■ ■ > Afc > . 

Give a partition A and a flag F,, the Schubert (sub)variety YxF, of G{k,n) is 

(1.1) YxF. := {E eG{k,n)\ dim En Fn^k+^^x,>^, i = I,-.., k}. 

This has codimension |A| := Ai -|- ■ • ■ -|- Afc in G{k, n). When A = □ := (1, 0, . . . , 0), 

YuF, = {E e G{k, n) \ diniF n > 1} , 

as the other conditions are redundant. We call □ a simple Schubert condition and Y]jF, 
a simple Schubert variety. It depends only upon Fn-k, so we also write Y\^Fn-k- All four 
Schubert varieties in the problem of four lines are simple. 

A Schubert problem is a list (A^, . . . , A™) of partitions with jA"*^! + ■ ■ ■ + |A™| = k{n—k). 
By Kleiman's Transversality Theorem [13], ii F^, . . . , F^ are general, then the intersection 

(1.2) YxiF^ n Yx2F^ n ■ ■ ■ n Yx-i^r 

is transverse and consists of finitely many /c-planes. The number d{X^ , . . . , X"^) of k- 
planes may be computed using the algorithms in the Schubert calculus (see [H] or [5] or 
the Introduction to [9]). The problem of four lines is an instance of the Schubert problem 
(□, □, □, □) in G(2, 4) and our analysis shows that d{D, □, □, □) = 2. 

We study Schubert problems in which all except possibly two Schubert conditions are 
simple. A simple Schubert problem on G{k, n) is one of the form 

(A, /i, □, □ ), 

k{n~k)~\\\~\^l\ 

where A, ix are not necessarily equal to □. We speak of the simple Schubert problem (A, /x) 
on G{k,n) (the k{n — k) — |A| — simple conditions □ are understood). 

The primary reason for limiting our study to simple Schubert problems in this paper 
is that these are Schubert problems that are complete intersections, and the off-the-shelf 
software that we use restricts us to complete intersections. 
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1.2. Galois groups of Schubert problems. According to Harris [8], Jordan [TT] showed 
how intrinsic structures of some enumerative problems could be understood in terms of 
Galois theory. Harris took the opposite approach — computing Galois groups of enumer- 
ative problems to expose the intrinsic structure of an enumerative problem. He showed 
that many enumerative problems have Galois group equal to the full symmetric group, 
demonstrating that these problems had no underlying structures. 

Harris' theory relating Galois groups to monodromy groups begins with a map 
V of degree d between irreducible complex algebraic varieties U and V. The function field 
C([/) of f/ is a degree d extension of the function field C{V) of V. These fields may be 
embedded into the field K of germs of meromorphic functions on a disc around a regular 
value V eV of f. If L G K is the normalization in K of the extension C{U)/C{V), then 
the Galois group G = GaA{L/C(y)) acts faithfully on the d points f~^{v), and this gives 
an embedding G "-^ Sd, where Sd is the symmetric group of the fiber f~^{y). 

Replacing U and V by Zariski open subsets if necessary, we may assume that the map 
f : U V \s a. degree d covering. A loop in V based at v has d hfts to f/, one for each 
point in the fiber f~^{v). Associating a point in the fiber f~^{v) to the endpoint of the 
corresponding lift gives a permutation in Sd- This defines the usual permutation action 
of the fundamental group of V on the fiber f~^{v). The monodromy group of the map 
f : U V is the image of the fundamental group of V in Sd- 

Proposition 1.1 (Harris [8]). For a map f : U ^ V as above, the monodromy group 
equals the Galois group. 

Given a Schubert problem (A^,...,A'") on G{k,n), let V be the space of m-tuples 
(Fj, . . . , F^) of fiags, a product of fiag manifolds. Define U to be the incidence variety 

(1.3) U ■-= {{H, . . . , F,™) G G{k, n)xV \ H e Y^^Fl for z = 1, . . . , m} . 

The fiber of U over a point H G G{k, n) is a product of the Schubert subvarieties 

{Fi\HeY^.Fi} 2 = l,...,m 

of the fiag manifold. Each of these is irreducible, and so U is irreducible. Let f : U ^ V 
be the other projection. Given v = [F^, . . . ,F^) G V, the fiber f~^{v) is the intersec- 
tion (11.21) . When v is general, this has d = d{X^, . . . , A"*) points, so that the map / has 
degree d- The Galois group of the Schubert problem is the Galois group of the extension 
C{U) /C{V). By Proposition 11.11 this is the monodromy group of the map f ■- U ^V- 

The point of this paper is that these monodromy groups may be computed using nu- 
merical homotopy continuation. For this, we first compute the points in a single fiber 
f~^{v)- Then, given a loop ip: [0,1] ^ V based at v {(p{0) = ip{l) = v), we numerically 
follow the points in the fibers f~^{ip{t)) as t runs from to 1. This computes the lifts of ip 
and thus the associated monodromy permutation. Computing sufficiently many of these 
monodromy permutations will enable us to recover the Galois group. While this gives the 
idea behind our method, we postpone more details until Section 13.11 
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2. HOMOTOPY CONTINUATION OF SIMPLE SCHUBERT PROBLEMS 

Homotopy continuation is a numerical method to compute all solutions to a system 
of polynomials given the solutions to a similar system. We use it to find all solutions 
to a simple Schubert problem and to compute elements of the monodromy group. We 
first describe the method of numerical homotopy continuation, then discuss polynomial 
formulations of Schubert problems, and finally explain the Fieri homotopy algorithm [21 
[TU] to find all solutions to simple Schubert problems. 

2.1. Homotopy continuation of polynomial systems. Suppose that we want to find 
all solutions to a 0-dimensional target system of polynomial equations 

(2.1) . . . = /2(Xi, . . . = ■■■ = f^iXi, . . . ,Xn) = 0, 

written as F{x) = 0. Numerical homotopy continuation finds these solutions if we have a 
homotopy, which is a system H{x,t) of polynomials in n+1 variables such that 

(1) The systems H{x, 1) = and F{x) = both have the same solutions; 

(2) We know all solutions to the start system H{x, 0) = 0; 

(3) The components of the variety defined by H{x,t) = include curves whose pro- 
jection to C (via the second coordinate t) is dominant; and 

(4) The solutions to the system H{x,t) = 0, where t G [0, 1), occur at smooth points 
of curves from (3) in the variety H{x, t) = 0. 

Given this, we restrict the variety H[x,t) = to t G [0,1] and obtain finitely many 
real arcs in x [0, 1] which connect (possibly singular) solutions of the target system 
H[x, 1) = to solutions of the start system H{x, 0) = 0. We then numerically trace each 
arc from t = to t = 1, obtaining all isolated solutions to the target system. 

The homotopy is optimal if every solution at t = is connected to a unique solution 
at t = 1 along an arc. This is illustrated in Figure [21 For simple Schubert problems, the 




H ^ H ^ H ^ 

10 10 1 

optimal not optimal not optimal 

Figure 2. Optimal and non-optimal homotopies 
Fieri homotopy algorithm is optimal. 

Remark 2.1. Homotopy continuation software often constructs a homotopy as follows. 
Let F(x) be the target system fl2.ll) and suppose we have solutions to a start system G{x). 
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Then for a number 7 G C with I7I = 1 define the linear homotopy 

H{x,t) := 7tF(x) + (1 -t)G(x). 

Then H{x, t) satisfies the definition of a homotopy for all but finitely many 7. The 
software detects the probability event that H{x, t) does not satisfy the definition when 
it encounters a singularity, and then it recreates the homotopy with a different number 7. 

Path following algorithms use predictor-corrector methods, which are conceptually sim- 
ple for square systems, where the number of equations equals the number of variables. 
Given a point t'-''-*) on an arc such that t'-''-* G [0, 1), the n x n matrix 



H, 



X 



is regular at {x^^\t^^^), which follows from the definition of the homotopy. Let Ht := 
{dHi/dt, dHjdtf. Given At, we set 

Ax := -Ati7,(xW,t(°))-ii7i(x(°),t(°)). 

For t(^) = + At, the point {x',t^^^) = + Ax,t^^^) is an approximation to the 
point (x^^\t^^^) on the same arc. This constitutes a first order predictor step. A cor- 
rector step uses the multivariate Newton method for the system H{x,t^^^) = 0, refining 
the approximate solution x' to a solution x^^\ In practice, the points x^^^ and x^^^ are 
numerical (approximate) solutions, and both the prediction and correction steps require 
that det 7^ at every point where the computation of the Jacobian matrix Hx is done. 

When the system is not square, additional strategies must be employed to enable the 
path following. Fortunately, simple Schubert problems are exactly the class of Schubert 
problems for which we have an optimal square homotopy (the Fieri homotopy). 

Cheater homotopies [20] are optimal homotopies constructed from families of polynomial 
systems. For example, given a Schubert problem (A^, . . . , A"^), let V be the space of all 
m-tuples (F,^, . . . , F^) of fiags. The total space of the Schubert problem 

U := {{H, . . . , F.'") e G{k, n)xV \ H e Y^Fl for z = 1, . . . , m} 

is defined by equations (see Section 12. 2p depending upon the point (F,^, . . . , F^^) G V. If 
ip: C ^ V is an embedding of C into V in which ip{0) and </?(!) are general m-tuples of 
flags and we write ^{t) = {F^{t), F^{t)), then 

= {{H,F^{t),...,F:rit))\H eY,.F:{t) for^ = l,...,m} . 

This is defined by a system H{x,t) = 0, which gives an optimal homotopy. We use 
this particular cheater homotopy to compute permutations in monodromy groups, called 
monodromy permutations. We give more details in Section 13. 1[ 

There, we describe the Fieri homotopy algorithm, which is a cheater homotopy where 
only one fiag in Lp{t) actually moves and the others remain fixed. The moving fiag is in 
general position when t = 1, but in a particular special position when t = 0, so that 
the Schubert problem becomes a union of other Schubert problems (whose solutions were 
previously computed and thus are known). 
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2.2. Equations for Schubert problems. Polynomial homotopy continuation methods 
require that our geometric problems are modeled by a system of polynomial equations. 
For efficiency, the number of variables should be minimized. We describe equations for 
Schubert varieties and then model Schubert problems by systems of equations which 
minimize the number of variables, stated in Proposition 12.21 below. 

Represent a fc-plane in C" as the row space of a by n matrix E with full rank and 
a flag by an invertible n by n matrix F, of constants, where the i-dimensional subspace 
in the flag is the row space of the flrst i rows of the matrix. The condition from (11.11) 
that dim_E fl Fn-k+i-Xi > i is 



(2.2) 



rank 



E 

Pn—k+i—\i 



< n — \i 



which is given by the vanishing of the determinants of all n+1— Aj by n+1— Aj submatrices 
of the n+i—\i by n matrix in fl2.2p . When Aj = the condition fl2.2l) is empty. 

Write S{E,F,,X) for the system consisting of these ) (n+r-A ) equations. 

The codimension |A| equals the number of equations only when A = □. In that case, 
S{E, F,, □) consists of the single equation 



(2.3) 



det 



E 



0. 



Since any two flags in general position are conjugate under a linear transformation, we 
always assume that two flags in (11.21) are flxed. Let the flag F, be deflned by setting Fi to 



be the span of e„, 



where ei. 



, Cr, form the standard basis for C". The 



Schubert variety YxF, has an open subset isomorphic to C*^'-" '^^ consisting of fc-planes 
that are the row space of an echelon matrix of the form 



0---0 1 * 

0---0 0---0 1 * 



* * • ■ ■ * 

* * • ■ ■ * 



0---0 0---0 0---0 1 

where 1 + Xk, 2 + Xk^i, . . . ,k + Xi are the columns with Is and * represents some number. 

To further reduce the number of variables, let the flag F^ be deflned by setting F- to 
be the span of ei, . . . , Cj. The skew Schubert variety [30] (or Richardson variety), 

Yx,^ := YxF.nV^F: 

has an open subset parameterized by matrices of the form 

1 * *0---0 

0---0 1 * *0 



E 



whose entries 



a j = i + Xk+i-i , 

if i + Xk+i-i < j <n 
otherwise . 



k + i — Hi 
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This parameterization is one-to-one if the product of the rightmost entries is non-zero, 

k 



i=l 



On the left below is E^a in G(2,4) and on the right is -E'210,110 in G{3, 7). 
[2A) 



1x00 
1 y 



1 a b c 
1 e 
1 /^ 



Given a Schubert problem A, /i, z/^, . . . , u"^, we will always take two of the flags to be 
these coordinate flags F, and F^, and consider intersections of the form 

Yx,^nY,iF^n---nY,..F:r. 

where the flags F^, . . . , F^ are general. By Kleiman's transversality theorem, all intersec- 
tions will lie in the subset of Yx^^ that is parameterized by matrices from Ex^^, and thus 
are solutions to the system of equations given by 

£iEx,„F^,i^'), 

This system is not necessarily square unless it is a simple Schubert problem. Since the 
homotopy continuation software we use is for square systems of polynomials, we restrict 
ourselves to simple Schubert problems. Write G for the n—k by n matrix Fn-k- Then 
S{Ex,fj,, F,, □) is a single equation (12. 3p that depends only on G. 

Proposition 2.2. A simple Schubert problem (A,yu) on G{k,n) is given by m := k{n — 
k) — \\\ — \ix\ matrices Gi, . . . , Gm each of size n—k by n, and the solutions are modeled 
by the system of equations 



(2.5) 



det 



E\^^ 
Gi 



det 



E 



Go 



det 



E 



X,li 



Gr^ 







(2.6) 



For example, the simple Schubert problem □, □ on G(2,4) is modeled by 



n n 1 

det 



det 



1x00 
1 y 

gu 912 913 9lA 
921 922 923 924 



9ii 9i2 9i3 



y 

9'lA 



where Gi = (gij) and G2 



921 922 923 924_ 

are matrices of constants. 



0, 



2.3. Pieri homotopy algorithm. We describe the simplified version of the Fieri ho- 
motopy algorithm [9], [10] that we use. The Fieri homotopy algorithm finds all solutions 
to those Schubert problems where all except possibly two partitions consist of a single 
part, (a, 0, . . . , 0). It is based on subtle geometric degenerations constructed in [27]. Both 
the algorithm and the degenerations enjoy a dramatic simplification for simple Schubert 
problems. The degenerations for these simple Schubert problems were introduced by 
Schubert 
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Example 2.3. Consider the simple Schubert problem (A, /i) = (210, 110) in G{3, 7). Begin 
with local coordinates (12. 4p for -^210,110 



E 



E. 



210,110 



1 a 6 c 
1 e 
1 /^ 



There are four columns, 1, 2, 3, and 6 not of the form n+i—k—^i. Let G be a general 
4-plane represented by a matrix in which these columns form a identity matrix 



(2.7) 



G 



1 

1 

1 





* * * 

* * * 

* * * 

* * 1 * 



Let G{t) be this matrix with each entry * scaled by t. Then 



(2.8) 



det 



E 
G{t) 



ceg + t(*c H h *cef) + f 



*bef)+t%*f + *af) 



where each * again represents a fixed number. 

When t = 0, the expression (12.81) becomes ceg. Let us investigate the consequences of 
ceg = 0. If we set c = in £'210,110, we get -E'210,210 and if we set g = 0, we get -^210,111, 



E' 



210,210 



1 a 6 
1 f/ e 
1 /^ 



E. 



210,111 



1 a 6 c 
1 e 
1 /0 



(If e = 0, then the row operation Ri ^ Ri — CR2 gives a matrix with c 
-E'210,210-) This computation in local coordinates shows that 



0, which lies in 



1^210,110 n YaG{0) 



Now suppose that we have general 4 by 7 matrices Gi, . 
the instance of the simple Schubert problem (210, 110): 



^210,210 U ^210,110 • 

,Gq,G7, and we wish to solve 



:2.9) 



det 



E 
Gi 



det 



E 

Ga 



det 



E 
Gj 







where G7 is the matrix (12. 7p . Replacing Gj by G{t) gives a homotopy of 7 equations in the 
coordinates a, . . . ,g where 6 equations are fixed (12. 9p and one depends on t (12. 8p . When 
t = 0, the latter becomes ceg = and the system splits into subsystems on -E'210,210 and 
-^210,111 involving the matrices Gi, . . . , Gq. Numerical continuation along this homotopy 
uses solutions to these smaller problems to obtain solutions to the system (12. 9p . 

Write A < if the components of the vector u — X are either or 1, with exactly one 1. 
For example, 110 < 210 and 110 < 111. Given partitions A,/i, define the (n— /c)-plane 



G, 



{ei\i^ {n-k+j-fij , for j = 1, . . . , k}} 
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Then if the non-zero entries of the matrix Ex^^ are aij, we have 



(2.10) 



det 



E 



X,fi 



G, 



n 



^i,n—k+i- 



which is the product of the rightmost non-zero entries in the rows of Ex,fj.- This determi- 
nant defines the charts Ex^,y for /i < z/. A child problem for the simple Schubert problem 
(A, fi) is one of the form (A, u) with /i < z/. 

The Fieri homotopy algorithm finds all solutions to a simple Schubert problem (A,/i) 
with fixed (but general) (n— A;)-planes Gi, . . . ,Gm where m -|- |A| + = k{n—k). We 
assume that we are given all solutions to all child problems (A, v) with fi < u and the 
(n— A;)-planes are Gi, . . . , Gm-i- If we let Gm(t) be a 1-parameter family of (n— A;)-planes 
with Gm(l) = Gm and G'm(O) = G^, then we obtain a homotopy 



det 



Ex,fi 
Gi 



det 



Ex,fi 

Gm-l 



det 



Ex,fi 

Gm{t) 



0, 



When t = this is the disjunction of child problems and when t = 1, it is the problem we 
wish to solve. 

This method recursively finds solutions to parent problems given solutions to their 
child problems. The depth of this recursion equals the dimension of the skew Schubert 
variety corresponding to the simple Schubert problem we wish to solve, i.e., the number 
of variables in the corresponding equations. The base case of this recursion is when 
|A| -|- = k{n—k), for then Ex^fj, is empty unless Aj + Hk+i-i = n—k for i = 1, . . . , fc, and 
in that case, Ex,^ gives the fc-plane spanned by {cj+A, \ i = 1, . . . ,k}. 

These homotopies are optimal. This is because they only follow solutions to the given 
Schubert problem and because the number d{\, /i) of solutions to a simple Schubert prob- 
lem satisfies the same recursion as the number of paths that are followed. Namely, if 
|A| -|- = k{n—k) then d{\, /i) = unless Aj -|- Hk+i-i = n—k ior i — 1, . . . , k, and then 
it equals 1. If |A| -|- < k{n—k), then 



d{X,fi) 



Remark 2.4. We do not quite use the homotopy we just described, as the equations 
involving t will in general have degree in t at least the minimum of k and n—k. Instead, 
we use the convex combination of the equations 



(2.11) 



7t det 



Ex,^l 

Gm 



t) det 



Ex,^l 
G„. 



0, 



Here, 7 G C has norm 1, I7I = 1. This has degree 1 in the homotopy parameter t. 

Doing this, the intermediate solutions will not necessarily be solutions to the Schubert 
problem, so we need to argue that this homotopy will remain optimal. For this, we appeal 
a little to the geometry of the Grassmannian. The equations we use define hyperplane 
sections of the Grassmannian in its Fliicker embedding, and the number of solutions d{\, fi) 
to the Schubert problem turns out to be the degree of the variety Ia,/^- Thus, replacing 
a family of hyperplanes defined by Schubert conditions (as in the Fieri homotopy) by 
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an arbitrary pencil of hyperplanes (12.111) will still give an optimal homotopy between 
solutions to the child problems and solutions to the parent problem. 

3. Description of software 

We divide the computation of the Galois group of a Schubert problem into three tasks. 

(1) Compute the solutions of a general instance of the problem, called a master set. 

(2) Use cheater homotopies to compute monodromy permutations of the master set. 

(3) Determine the group generated by the monodromy permutations. 

The first task may be accomplished by the brute-force application of a polynomial 
system solver. This is, however, inefficient. In Table [1], we compare the number of 
solutions to simple Schubert problems with A = yU = □ to the number of homotopy paths 
followed in polyhedral homotopies (as in the black-box solver of PHCpack). This is the 
volume of the associated Newton polytope and was computed with Polymake [7]. This 



k, n 


2,6 


2,7 


2,8 


2,9 


2,10 


2,11 


3,6 


3,7 


3,8 


4,6 


4,7 


5,7 


solutions 


14 


42 


132 


429 


1430 


4862 


42 


462 


6006 


14 


462 


42 


paths 


18 


67 


248 


919 


3426 


12843 


130 


3004 


74645 


42 


7156 


364 



Table 1 . Inefficiency of polyhedral homotopy for simple Schubert problems 

shows that the Fieri homotopy algorithm is an efficient alternative. 
Subsection 13.11 discusses the second task. 

The theory for the third task is beyond the scope of this paper. Based on preliminary 
computations, we conjectured that the Galois group is always the full symmetric group, 
and therefore we only check this. The fastest routine we have found to accomplish this is 
the isNaturalSymmetricGroup function of GAP [6]. 

In Section 13. 2^ we discuss implementations of our algorithm. Implementations and 
documentation of our computations are available at our website 



3.1. Computing monodromy permutations. Suppose that we have a master set of 
solutions to the simple Schubert problem (A, /i) on G{k, n), which is modeled by the system 
of equations 



(3.1) 



det 



det 



E 



Go 



det 



Gr^ 



0, 



where Gi, . . . ,Gm are fixed (n— /c)-planes. We follow these solutions along loops in the 
space of m-tuples of (n— /c)-planes to compute monodromy permutations. 

The off-the-shelf homotopy continuation software we use requires that the equations are 
linear in the homotopy parameter t, and so we follow piece-wise linear loops. For these, 
we fix all of the Gi except Gm, and exploit the linearity of the determinant in the row 
vectors gi, . . . , gn-k of Gm- If g'i is a vector not in Gm and we replace the row g^ by the 
convex combination (1 — t)gi + tg^, obtaining the pencil of planes Gm{t), then 



det 



E 



Gm{t) 



t) det 
t)Figi, 



Ex^fj, 



+ tdet 



E 



X,n 



Gm{l) 



tF{g 
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Replacing the equation in (13. ip involving Gm with this equation gives a linear homotopy 
between the system (13. ip and one with Grn(l) in place of Gm- 

Given a different {n—k)-p\a.iae G'^ spanned by g[, . . . , g'^_i^, we use these vectors to 
generate loops along which we can compute monodromy permutations. Suppose for il- 
lustration that n—k = 3 and the vectors are [a, b, c] for Gm and [a', b', c'] for G'j^. The 
different pencils that we may create correspond to the edges of a cube. 

[a',b',c'] 



[a', b', c] [a', b, d] [a, b', c'] 



[a', b, c] [a, b', c] [a, b, c'] 



[a, b, c] 

Our software offers 3 strategies to generate loops. 

• long loop goes from a vertex of the cube to its opposite and back, 

[a, b, c] — > [a', b, c] — > [a', b', c] — > [a, b', c'] 
— > [a, b', c'] [a, b, c'] — > [a, b, c] . 

• short loop uses only a square, 

[a, b, c] [a', b, c] [a' , b', c] — > [a, b' , c] — > [a, b, c] . 

• half loop makes use of just one edge, 

[a, 6, c] [a', 6, c] ^ [a, 6, c] , 

where the second homotopy is modified via a random number 7 G C, 

:= (1 -t)F(a',6,c) + -ftF{a,b,c). 

While we do not offer a proof that these loops will suffice to find all non-trivial permu- 
tations, we remark that they do suffice in the examples we considered. 

Example 3.1. Suppose that we have the simple Schubert problem (□, □) on G(2,4) as 
given by (12.61) with 



Gi 



G, 



-55 -8i 17+15i 40 + 99i -17-38i 

-67 + 25i -82 - 55i -99 - 80i -21 - 85i 

66 + 532 -73 - Ui 85 + 5i 67 + 16i " 

-53 - 85i 36 - 25i 2 + 81i -58 + 35i 



Its solutions nil and m2 are 

'1 -0.23714 - .0028980i 
1 



mi 





.51680 - .10520Z 



m2 



1 .97009 + 1.2705Z 

1 .44336 + .38248^ 
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The short loop strategy with 

, _ r 33 - 84i 21 -i 59 + 9Ai -94 + 89i " 
[ -15 - 19i 29i 79 + 51i 89 + 3i 

creates a non-trivial monodromy permutation. The paths followed during the homotopies 
are drawn in Figure [31 where the large circles are the values taken at the endpoints of the 
different homotopies. 




Second coordinate in position 2, 4 



Figure 3. Paths tracked in Example [XT! 

3.2. Implementation. We have two Maple implementations of our algorithms using the 
package PHCmaple [17j to interface with PHCpack |33j, which performs the numerical 
polynomial homotopy continuation. PHCmaple produced the graphic of Figure [31 The 
second implementation may alternatively call Bertini [1]. 

Our prototype implementation was carried out entirely in Maple to take advantage of 
Maple packages to generate the equations and to manage the monodromy group, while 
using the black-box solver in PHCpack to compute the master sets of solutions. The largest 
problem this implementation could treat was the simple Schubert problem (210, 200) on 
G(3, 7) — it showed that the Galois group is the full symmetric group 5*91. Previously, the 
largest Schubert problem whose Galois group that was proven to be the full symmetric 
group was the simple Schubert problem (20, 10) on G(2, 6) with 9 solutions [3]. 

Our second implementation also uses Maple and either PHCpack or Bertini. However, 
it relies on the Fieri homotopy algorithm to compute master sets of solutions and GAP 
to manage the monodromy groups, removing the two main computational bottlenecks 
of the prototype. The largest problems this implementation has treated are (2100, □) 
in G{4:, 8) with 8580 solutions, (210, 210) in ^(3, 8) with 10329 solutions and (210, 200) 
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on G(3, 9) with 17589 solutions. It showed that all of these have Galois group equal to 
the full symmetric group. The computation for the problem (□, □) on G(3, 8) with 6006 
solutions is the basis of the Numerical Theorem. 

We computed Galois groups of the simple Schubert problems (□, □) on all small Grass- 
mannians, using the short loops strategy. They were run on several different computers, 
including an AMD Athlon 64 Dual Core Processor 4600+ with CPU clock speed of 2400 
MHz and 1 GB of memory whose timings (using PHCpack) are reported in Tabled These 



k, n 


2,4 


2,5 


2,6 


2,7 


2,8 


2,9 


2,10 


solutions 


2 


5 


14 


42 


132 


429 


1430 


time 


12s 


27s 


19s 


51s 


4.2m 


20.5m 


2.6h 


permutations 


4 


6 


5 


6 


7 


4 


7 



k, n 


3,5 


3,6 


3,7 


3,8 


3,9 


4,6 


4,7 


4,8 


solutions 


5 


42 


462 


6006 


17589 


14 


462 


8580 


time 


12s 


35s 


17.9m 


18.6h 


78.2h 


15s 


23.5m 


44.5h 


permutations 


4 


4 


5 


6 


7 


5 


5 


7 



seconas, m : 



minutes, h := hours 
Table 2. Timings of Galois group computation 



reported times are not CPU times, but actual elapsed (wall clock) time, and so may ex- 
ceed CPU time by 10 to 20 %. We also record the number of permutations we needed to 
compute. The entry in G(3, 9) is the Schubert problem (210, 200) and the entry in G*(4, 8) 
is the Schubert problem (2100, □). 

We ran some of these computations in Bertini. It was unable to compute examples in 
more than 10 variables, and was markedly slower for the largest computations it com- 
pleted, on G(2,8), G(3, 7), and (7(4,7). On the other hand, Bertini provided an inde- 
pendent verification that the Galois groups were indeed the full symmetric groups. We 
remark that Bertini is new software and its efficiency will likely improve. 

4. Conclusions and Future work 

This paper demonstrates the feasibility of homotopy continuation as tool to study the 
Galois groups of enumerative problems. It also implicitly provides several challenges to 
the numerical homotopy community. Perhaps the most serious is the current lack of 
certifiability of computations in numerical homotopy software. While numerical methods 
will increasingly outperform symbolic algorithms in algebraic geometry, they are currently 
inferior in that their results do not come with certification. Certificates for numerical com- 
putations do exist in theory, for example in Shub and Smale's [23] alpha-theory, and there 
is a need for their implementation. In fact, even more reliable numerical techniques, such 
as the interval step control proposed in [12], are sidestepped by homotopy continuation 
software developers mostly due to the perceived complexity of implementation and the 
expected slower performance in comparison with heuristic methods. Robust, off-the-shelf 
software to handle polynomial systems that are not complete intersections is also needed 
to deal with ideals in algebraic geometry, which are typically not complete intersections. 
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Two further theoretical problems are not addressed in this paper. While the sampling 
of the fundamental group of the base space provided by the loop-generating heuristics of 
subsection 13.11 is sufficient to generate the Galois group in the examples we considered, 
a better understanding of the topology of a complement of an algebraic variety in a 
Grassmannian is needed to prove that this sampling is always sufficient. Second, we do 
not know how to certify that the computed set of permutations generates the whole Galois 
group, if it is not the full symmetric group. 

The tools that we use could be applied more systematically to other problems in enu- 
merative geometry. To this end, we plan a comprehensive project exploring the limit of 
computability of Galois groups of Schubert problems along several fronts. This will involve 
the software and algorithms described here — perhaps also incorporating H0M4PS [15]. 
We will also write parallel software implementing algorithms to compute Galois groups 
of Schubert problems that are not compete intersections, as well as pushing the limits of 
the symbolic methods of Billey and Vakil [3] and of Vakil's combinatorial algorithm [52] . 
This will be a large and distributed computation as in [21| which should give a catalog of 
the Galois groups of several tens of thousands of Schubert problems. 
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